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Abstract 

This article deals with the numerical resolution of Markovian backward stochastic differential 
equations (BSDEs) with drivers of quadratic growth with respect to z and bounded terminal con- 
ditions. We first show some bound estimates on the process Z and we specify the Zhang's path 
regularity theorem. Then we give a new time discretization scheme with a non uniform time net for 
such BSDEs and we obtain an explicit convergence rate for this scheme. 

1 Introduction 

Since the early nineties, there has been an increasing interest for backward stochastic differential equa- 
tions (BSDEs for short). These equations have a wide range of applications in stochastic control, in 
finance or in partial differential equation theory. A particular class of BSDE is studied since few years: 
BSDEs with drivers of quadratic growth with respect to the variable z. This class arises, for example, in 
the context of utility optimization problems with exponential utility functions, or alternatively in ques- 
tions related to risk minimization for the entropic risk measure (see e.g. |[T3l ). Many papers deal with 
existence and uniqueness of solution for such BSDEs: we refer the reader to ifTTlfTSl when the terminal 
condition is bounded and EllUlQl for the unbounded case. Our concern is rather related to the simulation 
of BSDEs and more precisely time discretization of BSDEs coupled with a forward stochastic differen- 
tial equation (SDE for short). Actually, the design of efficient algorithms which are able to solve BSDEs 
in any reasonable dimension has been intensively studied since the first work of Chevance ||6l, see for 
instance ifTOl fTl fTTl . But in all these works, the driver of the BSDE is a Lipschitz function with respect to 
z and this assumption plays a key role in theirs proofs. In a recent paper, Cheridito and Stadje [iS] stud- 
ied approximation of BSDEs by backward stochastic difference equations which are based on random 
walks instead of Brownian motions. They obtain a convergence result when the driver has a subquadratic 
growth with respect to z and they give an example where this approximation does not converge when the 
driver has a quadratic growth. To the best of our knowledge, the only work where the time approxima- 
tion of a BSDE with a quadratic growth with respect to z is studied is the one of Imkeller and Reis |[T4l . 
Let notice that, when the driver has a specific forrrQ, it is possible to get around the problem by using 
an exponential transformation method (see lfT5l ) or by using results on fully coupled forward-backward 
differential equations (see Q). 



'Roughly speaking, the driver is a sum of a quadratic term z i— >■ C\z\^ and a function that has a linear growth with respect 



to 2. 



1 



] INTRODUCTION 



2 



To explain ideas of this paper, let us introduce {X, Y, Z) the solution to the forward backward system 

Xt = x + [ b{s,Xs)ds+ [ a{s)dWs, 
Jo Jo 

Yt = 9{Xt) + J^ f{s,Xs,Ys,Zs)ds- ZsdWs, 

where g is bounded, / is locally Lipschitz and has a quadratic growth with respect to z. A well-known 
result is that when 5 is a Lipschitz function with Lipschitz constant Kg, then the process Z is bounded 
by C{Kg + 1) (see Theorem 13. lb . So, in this case, the driver of the BSDE is a Lipschitz function with 
respect to z. Thereby, a simple idea is to do an approximation of {Y, Z) by the solution (F^, Z'^) to the 
BSDE 

Yi" = gM{XT) + 1^"^ f{s,Xs,Yj',Z^)ds- Z^dWs, 

where is a Lipschitz approximation of g. Thanks to bounded mean oscillation martingale (BMO 
martingale in the sequel) tools, we have an error estimate for this approximation: see e.g. |[T4l |2l or 

— a 

Proposition 14.21 For example, if g is a-Holder, we are able to obtain the error bound CKg~" (see Propo- 
sition 14.101 ). Moreover, we can have an error estimate for the time discretization of the approximated 
BSDE thanks to any numerical scheme for BSDEs with Lipschitz driver. But, this error estimate de- 
pends on Kgj^: roughly speaking, this error is Ce 3Nn~ with n the number of discretization times. 
The exponential term results from the use of Gronwall's inequality. Finally, when g is a-Holder and 
Kg^ = N, the global enw bound is 

Ci^ + - . (1.1) 

yN~ n J 

So, when N increases, will have to become small very quickly and the speed of convergence turns 
out to be bad: if we take N = logn)^^^ with < e < 1, then the global enw bound becomes 

Cs (log n) . The same drawback appeai^s in the work of Imkeller and Reis HH. Indeed, their idea 
is to do an approximation of (Y, Z) by the solution {Y^ , Z^) to the truncated BSDE 

Y,^ = g{XT) + f{s,Xs,Y^'',hN{Z^))ds- j\^dWs^ 

where : M^^*^ — )• M^^*^ is a smooth modification of the projection on the open Euclidean ball of 
radius N about 0. Thanks to several statements concerning the path regularity and stochastic smoothness 
of the solution processes, the authors show that for any /? ^ 1, the approximation error is lower than 
CjjN~l^ . So, they obtain the global error bound 

and, consequently, the speed of convergence also turns out to be bad: if we take N = log with 
< e < 1, then the global error bound becomes C/j^e (logn)"''^^. 

Another idea is to use an estimate of Z that does not depends on Kg. So, we extend a result of ISl 
which shows 

l^*l^^i+ (r!t)V2 ' 0^*<^- (1-3) 
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Let us notice that this type of estimation is well known in the case of drivers with linear growth as a 
consequence of the Bismut-El worthy formula: see e.g. lITOl . But in our case, we do not need to suppose 
that a is invertible. Then, thanks to this estimation, we know that, when t < T, /(i, ., ., .) is a Lipschitz 
function with respect to z and the Lipschitz constant depends on t. So we are able to modify the classical 
uniform time net to obtain a convergence speed for a modified time discretization scheme for our BSDE: 
the idea is to put more discretization points near the final time T than near 0. The same idea is used 
by Gobet and Makhlouf in fT2l for BSDEs with drivers of linear growth and a terminal function g not 
Lipschitz. But due to technical reasons we need to apply this modified time discretization scheme to the 
approximated BSDE: 

Jt Jt 

with 

f%s, X, y, z) := ls<T~ef{s, X, y, z) + ls^T~ef{s, X, y, 0). 

Thanks to the estimate (II. 3I ). we obtain a speed convergence for the time discretization scheme of this 
approximated BSDE (see Theorem 14. 8 1 ). Moreover, BMO tools give us again an estimate of the approx- 
imation error (see Proposition 14.21) . Finally, if we suppose that g is a-Holder, we prove that we can 
choose properly N and e to obtain the global error estimate Cn (2-a)(2+K)~2+2a (^gg Theorem 14.131) 
where K > depends on constant M2 defined in equation (11.31 ) and constants related to /. Let us notice 
that such a speed of convergence where constants related to /, g, b and a appear in the power of n is un- 
usual. Even if we have an eiTor far better than (11.11 ) or (11.21 ). this result is not very interesting in practice 
because the speed of convergence strongly depends on K. But, when b is bounded, we prove that we can 
take M2 as small as we want in (11.31) . Finally, we obtain a global error estimate lower than C^n~*^"~^\ 
for all 77 > (see Theorem l4. 161 ). 

The paper is organized as follows. In the introductory Section 2 we recall some of the well known 
results concerning SDEs and BSDEs. In Section 3 we establish some estimates concerning the process 
Z: we show a first uniform bound for Z, then a time dependent bound and finally we specify the classical 
path regularity theorem. In Section 4 we define a modified time discretization scheme for BSDEs with a 
non uniform time net and we obtain an explicit eiTor bound. 

2 Preliminaries 
2.1 Notations 

Thi^oughout this paper, (Wf)t>o will denote a d-dimensional Brownian motion, defined on a probability 
space {U,T,F). For t ^ 0, let J^t denote the cr-algebra a{Ws; ^ s ^ t), augmented with the P-null 
sets of J^. The EucUdian norm on M'^ will be denoted by | . | . The operator norm induced by | . | on the 
space of Unear operator is also denoted by | . | . For p ^ 2, m G N, we denote further 

• 5^(M™), or when no confusion is possible, the space of all adapted processes (yt)fg[o^r] with 
values in M^normed by ||y||5p = E[(supig[o_r] 1^*1)^]^/^; 5°° (K'"), or 5°°, the space of bounded 
measurable processes; 

• A^^(M"^), or Ai^, the space of all progressively measurable processes (^t)jg[o,T] with values in 

R"" normed by \\Z\\j^p = E[{J^ \Zsf ds)P/2]i/p. 

In the following, we keep the same notation C for all finite, nonnegative constants that appear in our 
computations: they may depend on known parameters deriving from assumptions and on T, but not on 
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any of the approximation and discretization parameters. In the same spirit, we keep the same notation t] 
for all finite, positive constants that we can take as small as we want independently of the approximation 
and discretization parameters. 

2.2 Some results on BMO martingales 

In our work, the space of BMO martingales play a key role for the a priori estimates needed in our 
analysis of BSDEs. We refer the reader to |[T6l for the theory of BMO martingales and we just recall 
the properties that we will use in the sequel. Let = /q (psdWg, t G [0, T] be a real square integrable 
martingale with respect to the Brownian filtration. Then $ is a BMO martingale if 



\HbMO= sup E[($)t- (^)rl-^r]'/' = SUp E 
rg[0,T] Te[0,T] 



T 



1/2 



< +00, 



where the supremum is taken over all stopping times in [0, T]; (<I>) denotes the quadratic variation of 
In our case, the very important feature of BMO martingales is the following lemma: 

Lemma 2.1. Let ^ be a BMO martingale. Then we have: 

L The stochastic exponential 

£{^)t = St = exp (^j^ <i)sdWs - ^ \(l)sf ds^ , ^ t ^ T, 

is a uniformly integrable martingale. 

2. Thanks to the reverse Holder inequality, there exists p > 1 such that £t G L'^- The maximal p with 
this property can be expressed in terms of the BMO norm of^. 



3. Mn G N*, E 



2.3 The backward-forward system 

Given functions b, a, g and /, for x G M'^ we will deal with the solution {X, Y, Z) to the following 
system of (decoupled) backward-forward stochastic differential equations: for t G [0, T], 

Xt = x + [ b{s,Xs)ds+ [ a{s)dWs, (2.1) 
Jo Jo 

Yt = g{XT)+ [ f(s,Xs,Ys,Zs)ds- [ Z^dW^. (2.2) 
Jt Jt 

For the functions that appear in the above system of equations we give some general assumptions. 

(HXO). b : [0, T] X M'^ R'^, a : [0, T] R'^'^'^ are measurable functions. There exist four positive 
constants Mb, Kb, M^ and such that yt, t' G [0, T], \/x, x' G R'^, 

\b{t,x)\ ^ Mb{l + \x\), 

\b{t,x)-b{t',x')\ ^ Kb{\x-x'\ + \t-t'\^''^), 

^ M,, 

\a{t)-a{t')\ ^ K^\t-t'\. 
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(HYO). / : [0,r] X M'^ X M X Ri^'^ M, 5 : M'' M are measurable functions. There exist 
five positive constants M/, Kf^^, Kf^y, Kj^^ and Mg such that Vt G [0,r], \/x,x' £ W'-, \/y,y' G M, 

\f{t,x,y,z)\ ^ M;(l + |y| + |z|2), 
- f{t,x\y\z')\ ^ A'/^^, |x - x'\ + A/^y |y - y'\ + {Kf^^ + Lf^^{\z\ + \z'\)) \z - z'\ , 
\g{x)\ ^ Mg. 

We next recall some results on BSDEs with quadratic growth. For their original version and their proof 
we refer to ifTTIl. Ell and 

Theorem 2.2. Under (HXO), (HYO), the system f[Zll-f[Z2\l has a unique solution (X, Y, Z) € S'^xS°°x 
A4'^. The martingale Z *W belongs to the space ofBMO martingales and \\Z * VF||bmo '^^^y d^P^^ds 
on T, Mg and Mf. Moreover, there exists r > 1 such that £{Z * W) G L^. 

3 Some useful estimates of Z 
3.1 A first bound for Z 

Theorem 3.1. Suppose that (HXO), (HYO) hold and that g is Lipschitz with Lipschitz constant Kg. Then, 
there exists a version of Z such that, Vt G [0, T], 

\Zt\ ^ e^^^^+^i'y^'^M,{Kg + TKf^,). 

Proof. Firstly, we suppose that h, g and / are differentiable with respect to x, y and z. Then {X, Y, Z) 
is differentiable with respect to x and [VX, VY, VZ) is solution of 

VXt = Id + [ Vb{s,Xs)VXsds, (3.1) 
Jo 

VYt = VgiXT)VXT-J^ VZsdWs (3.2) 

+ VJis, X,,Y,, Zs)VX, + Vyfis, X,,Y,, Zs)VY, + V,/(s, X„ Z,)VZ,ds, 



where VXt = {dXl/dx^)i^,j^d, VF* = \dYt/dx^)i^j^d e R^""^, V Zt = {dZi/dx^)i^ij^d and 
VZsdWs means 



[ (vz.ydwi 



with {SI Zy denoting the i-th line of the d x d matrix process V Z. Thanks to usual transformations on 
the BSDE we obtain 

g/o ^yf{s,Xs,Ys,Zs)ds^^Y^ = ^v-f^''^^='^"^''')'^''\/g(XT)VXT - / 6-^0 Vy/(«,X„,Y„,Z„)d«y 

rT 

+ ei'o^«/("'^-^".^")'^"V,./(s,X„n,Z,)VX,ds, 
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with dWs = dWs - VJ{s, X^, Ys, Zs)ds. We have 

2 



VJ{s,Xs,Ys,Zs)dWs 







BMO 



sup E 

r6[0,T] 



\V,f{s,Xs,Ys,Zs)\'ds 











( 1 + sup E 


\ \Zsi\ ds 






V ^e[o,r] 









C 1 + IIZ* VFI 



BMO 



Since Z *W belongs to the space of BMO martingales, || ^/^/(s, X^, Yg, 



^Wbmo 



< +00. 



Lemma IZT] gives us that £{Jq Vzf{s, Xs,Ys, Zs)dWs)t is a uniformly integrable martingale, so we ar^e 
able to apply Girsanov's theorem: there exists a probability Q under which (W^)tg[o,T] is a Brownian 
motion. Then, 



J,lS7yfis,Xs,Ys,Zs)ds^Y^ 



E^ 



Jo ^yf('''^='Y='^=)'i''Vg{XT)VXT 



and 



ivy*! ^ e(^^^+^/-'^)^(i^3 + TKf^,), (3.3) 

because | VX^I ^ e^^'^ . Moreover, thanks to the Malliavin calculus, it is classical to show that a version 

of {Zt)t(i[Q^T] is given by (Vyt(VXt)- V(t))ig[o,T]- So we obtain 

\Zt\ ^ e^^'^M, \VYt\ ^ e^''^^+^i'y'^^M,{Kg + TKj^,), a.s., 

because |VXj"^| ^ e^^T^ 

When b, g and / are not differentiable, we can also prove the result by a standard approximation and 
stability results for BSDEs with Unear growth. □ 

3.2 A time dependent estimate of Z 

We will introduce two alternative assumptions. 

(HXl). b is differentiable with respect to x and a is differentiable with respect to t. There exists 
A G M+ such that Vr/ G R'^ 



*7?a(s)[*a(s)*V6(s,x) - A |*r/^7(s)|^ 



(3.4) 



(HXl'). 0- is invertible and Vt G [0, T], \(j{t)~^\ ^ M^-i. 



Example. Assumption (HXl) is verified when, Vs G [0, T], Vb{s, .) commutes with a{s) and 3^ : 
[0, T] R'^'"^ bounded such that a'{t) = a{t)A{t). 

Theorem 3.2. Suppose that (HXO), (HYO) hold and that (HXl) or (HXl ') holds. Moreover, suppose that 
g is lower ( or upper) semi-continuous. Then there exists a version of Z and there exist two constants 
C, C G M"*" that depend only in T, Mg, Mf, Kf .j., Kf y, Kf^^ and Lf^^ such that, Vt G [0, T[, 

\Zt\ i^C + C'{T -t)-^/^. 
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Proof. In a first time, we will suppose that (HXl) holds and that f,gaie differentiable with respect to 
X, y and z. Then {Y, Z) is differentiable with respect to x and {VY, VZ) is the solution of the BSDE 

VYt = Vg{XT)VXT- j\zsdWs 

+ ^xfis, Xs,Ys, Zs)VXs + Vyfis, Xs, Ys, Zs)VYs + VJ{s, X^, Y,, Zs)VZsds. 
Thanks to usual transformations we obtain 

^0 

e/o^v,/(s,x.,n,z,)dsy^^^^^y^^ ^ f e-^o^ynn,Xu,Yu,z^)duy^j(^^^j^^^Y„Zs)VXsds 

Jo 

f-T 

_ j e^O '^yf(u,Xu,Yu,Zu)du^2^yy 



with dWs = dWs — Vzf{s,Xs,Ys, Zs)ds. We can rewrite it as 



Ft = FT 



(3.5) 



with 



:= e^o ^yf{s,XsXs,Zs)dsyYt + / e^o v,/(«,x„,y„,z„)d«y^y.(g^ js^:^^ y^^ Z,)VX,(is. 

JO 

Z*W belongs to the space of BMO martingales so we are able to apply Girsanov's theorem: there exists 
a probability Q under which (l^)tg[o,T] is a Brownian motion. Thanks to the Malliavin calculus, it is 
possible to show that (Vlf(VXf)~^(T(i))(g[o^r] is a version of Z. Now we define: 



at 



f e/o ^yf(u,x^X^,z^)duyj^^^ x^^y^^ Z,)VXsds{VXtr^a{t), 
Jo 



Ft ■■= e^'Ft(VXt)~\ 

Since dVXt = Vb{t, Xt)VXtdt, then (i(VXj)~i = -(VXt)~iV6(i, Xt)(it and thanks to Ito's for- 
mula, 

dZt = dFt{VXt)-^a{t) - Ft{VXt)~^Vh{t,Xt)a{t)dt + Ft{V Xt)'^ a' {t)dt , 

and 

d{e^^Zt) = Ft{\Id - Vb{t,Xt))a{t)dt + Fta'{t)dt + e^*dFj(VXj)-V(t). 

Finally, 



A 



Fta{t) - Fta{t)['a{tYVb{t,Xt) - 'a'{t)fFt 



dt + dM^, 



with Mt := /q e^'*dFs(VXs)"V(s) and A// a Q-martingale. Thanks to the assumption (HXl) we are 

2 



able to conclude that 



e^'Zt 



r 




2 








ds 













is a Q-submartingale. Hence, 

2 



> e 



2At 



Zt 



(T-t) 



> e 



2At 



Jo'^yf{s,Xs,Ys,Zs)ds^ _^ 



t + a* 



(T — t) a.s., 
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which impUes 

\Zt\^iT-t) 



-2Xt -2 J^Vyf{s,X,,Ys,Zs)ds 2\t 



Ji^yfis,Xs,Ys,Zs)ds^ + 



at - at 



(T-t) 



^ C e 



2Xt 



J^\7yf{s,Xs,Ys,Zs)ds2^ + 



at 



+ 1 (T-t) 



^ C E' 



r 




2 






Zs 


ds 


Ft 











(T-t) 



a.s. 



with C a constant that only depends on T, K^, Kf^y and A. Moreover, we have, a.s., 

rT 



r r'^ 




2 










Zs 


ds 


Ft 


^ CE^ 

















iZoP + la.l^ ds 



Ft 



^ C \\Z 



\BMO( 



+ iT-t)). 

But ||^||ba/o(iQ) '^'^^^ "^^^ depend on Kg because {Y, Z) is a solution of the following quadratic BSDE: 

Yt = g{XT) + r {f{s, Xs, Ys, Zs) - ZsVJis, X„ F,, Zs)) ds - ZsdWs- (3.6) 



Jt Jt 

Finally ^ C (l + (T - t)-^'"^) a.s.. 

When a is invertible, the inequality (13.41 ) is verified with A := M^^i{M^Kb + Ka)- Since this 
A does not depend on V6 and a', we can prove the result when b{t, .) and a are not differentiable by a 
standard approximation and stability results for BSDEs with linear growth. So, we are allowed to replace 
assumption (HXl) by (HXl'). 

When / is not differentiable and g is only Lipschitz we can prove the result by a standard approxi- 
mation and stability results for linear BSDEs. But we notice that our estimation on Z does not depend on 
Kg. This allows us to weaken the hypothesis on g further: when g is only lower or upper semi-continuous 
the result stays true. The proof is the same as the proof of Proposition 4.3 in ||8l. □ 

Remark 3.3. The previous proof gives us a more precise estimation for a version of Z when f is differ- 
entiable with respect to z: \/t G [0, T], 



\Zt\ (^C + C'E^ 



\ZJ'^ ds 



Ft 



1/2 



(T-t) 



-1/2 



Remark 3.4. When assumptions (HXl ) or (HXl ') are not verified, the process Z may blow up before T. 
Zhang gives an example of such a phenomenon in dimension 1: we refer the reader to example 1 in A20I/ . 

3.3 Zhang's path regularity Theorem 

Let = to < *i < • • • < = ^ be any given partition of [0, T], and denote (5„ the mesh size of this 
partition. We define a set of random variables 



1 



-E 



Z.,ds 



U 



ti+l — ti 

Then we are able to precise Theorem 3.4.3 in 1(211 : 



Ft. 



Vie{0, ...,n-l}. 



Theorem 3.5. Suppose that (HXO), (HYO) hold and g is a Lipschitz. function, with Lipschitz constant Kg. 
Then we have 



n-1 
1=0 



Zt-Zt. 



dt 



where C is a positive constant independent of 6n and K„ 
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Proof. We will follow the proof of Theorem 5.6., in |[T4l : we just need to specify how the estimate 
depends on Kg. Firstly, it is not difficult to show that Zt- is the best 7"^- -measurable approximation of Z 



E 



ti + l 



Zt - Zt^ dt 



inf E 



ti+i 



Zt — Zi\ dt 



In particular. 



E 



ti+i 



Zt - Zt, dt 



^ E 



\Zt - Zti\ dt 



In the same spirit as previous proofs, we suppose in a first time that h, g and / are differentiable with 
respect to x, y and z. So, 

Zt - Zt^ = VYt{VXtr^a{t) - VYt^{VXt,)-^a{ti) = h + h + h, a.s., 

with h = VYt{VXt)-\ait) - a{ti)), h = yYt{{VXt)~^ - (VXtJ-i)c7(ti) and h = "^{Yt - 
Yt^){y Xt^)^^ a{ti) . Firstly, thanks to the estimation (13.3b we have 

^ \VYt\' e^^^'^Kl - ^ C{1 + Kl)5l. 

We obtain the same estimation for I/2I because 

\{vxtr^ - (vxtj-^i ^ [\vxsr'vb{s,Xs)d 

Jti 

Lastly, I/3K M.e^"^ \VYt - Vr^J. So, 

n-l 



Ki,e^^'^ \t-t 



j=0 





n.-l 


r^'\h\^dt 


^ C5n^E 




i=0 



ess sup \VYt-VYt^\ 

t(^[ti,ti + l] 



By using the BSDE (1X21) . (HYO), the estimate on VXg and the estimate (1331) . we have 

\VYt-VYtf 

^c(^l^^{Cil + Kg) + \VJ{s,Xs,Y,,Zs)\\VZs\)ds^ +c(^j^ VZsdw)j . 



The inequalities of Holder and Burkholder-Davis-Gundy give us 

n-l 

ess sup \VYt-VYt^\ 

1 / „+. 



j=0 



^ C(1 + k2) + C J^E 

i- 

^ C{\ + K]) + C7E 



|V,/(s,X„n,Z,)||VZ,|(isj +CE 

^ (l + |Z,|2)(is^ (^^ \VZ,\^ds^+j^ \VZs\^ds 

\VZ.,\^ds 



\VZs\ ds 



\ZJ ds 



4 CONVERGENCE OF A MODIFIED TIME DISCRETIZATION SCHEME FOR THE BSDE 1 



for all p > 1 and g > 1 such that 1/p + 1/q = I. But, {VY, VZ) is solution of BSDE ^2^, so, from 
Corollary 9 in [2], there exists q that only depends on \\Z * W\\^j^,jq such that 



E 



T 



|VZ,| ds 



1/9 



^ C{1 + K\ 



Moreover, we can apply Lemma |2?T] to obtain the estimate 



E 



T 



\ZJ\ ds 



i/p 



Finally, 



and 



n-l 
i=0 



n-l 



\Zt-Zt}\ dt 



\hrdt 



^ C(l + K^M 



n—1 r f'f- 
1=0 

^ C{l + Kl)6n. 



h\^ + \l2\^ + \h\^]dt 



□ 



4 Convergence of a modified time discretization scheme for the BSDE 
4.1 An approximation of the quadratic BSDE 

In a first time we will approximate our quadratic BSDE (12.21 ) by another one. We set e s]0, T[ and 
N en. Let {Yf'", Z^'") the solution of the BSDE 



with 



Yf'' = g^iXr) + r fis, X„ Z^'')ds - f Z^'^dWs 
Jt Jt 



fis, X, y, z) := lsi:T-efis, x, y, z) + ls>T-ef{s, X, y, 0), 



(4.1) 



and 5 AT a Lipschitz approximation of g with Lipschitz constant N. verifies assumption (HYO) with the 
same constants as /. Since gj\[ is a Lipschitz function, Z^'^ has a bounded version and the BSDE (14.11) 
is a BSDE with a linear growth. Moreover, we can apply Theorem l3.2l to obtain: 



Proposition 4.1. Let us assume that (HXO), (HYO) and (HXl) or (HXl ') hold. There exists a version of 
Z^'^ and there exist three constants Mz^i, Mz^2, Mz^s G that do not depend on N and e such that, 
VsG[0,T], 



M. 



z,2 



A(M,,3(A^+1)). 



(T- s)i/2_ 

Thanks to BMO tools we have a stabihty result for quadratic BSDEs (see Q and HH): 
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Proposition 4.2. Let us assume that (HXO) and (HYO) hold. There exists a constant C that does not 
depend on N and e such that 



E 



sup 

te[o,T] 



Y^'^ - Y, 


2 


+ E 


if 


ry-N,e ry 


2 

dt 








Jo 







with 



ei(iV):=E \g^{XT) - 9{Xt , 



|2g 



^C(ei(iV) + e2(A^,e)) 

1/9 



,1 



e2(iV,e) :=E 
and q defined in Theorem^ 



2q 



1/9 



Then, in a second time, we will approximate our modified backward-forward system by a discrete- 
time one. We will slightly modify the classical discretization by using a non equidistant net with 2n + 1 
discretization times. We define the n + 1 first discretization times on [0, T — e] by 



tk = T[l 



e \ fc/" 

t) 

and we use an equidistant net on [T — e, T] for the last n discretization times: 

'2n - 



tk = T 



n 



e, n ^ k ^ 2n. 



We denote the time step by {hj. := t^j^i — ife)o^fe^2n-i- We consider (^t" )o^fc!j2n the classical Euler 
scheme for X given by 



X 



We denote ps 



XI + hub{tu,Xl) + a{tu){Wt,^, -Wt,), < A: < 2n - 1. 
the projection on the ball 

B ( 0, 1 



(4.2) 



2,2 



(r-s)V2 



with M^ i and M2 2 given by Proposition 14.11 Finally we denote (y^.^,"^ ^^-e.'^) our time approx- 
imation of This couple is obtained by a slight modification of the classical dynamic 
programming equation: 



N.e.n 



Y, 



Ptu+. (ylft, Kf;r(^*.H-i - J]) ' ^ A; ^ 2n - 1, (4.3) 
E*. Kf;r ] + [f{t,, XI , y,f;f , , ^ ^ 2n - 1, (4.4) 



■^The authors of 1141 obtain this result with instead of q. Nevertheless, we are able to obtain the good result by applying 
the estimates of (2). 
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where E^j. stands for the conditional expectation given Tt^- Let us notice that the classical dynamic 
programming equation do not use a projection in (14.31 ): it is the only difference with our time approx- 
imation, see e.g. ifTTTl for the classical case. This projection comes directly from the estimate of Z in 
Proposition 14.11 The aim of our work is to study the error of discretization 



e{N,e,n) := sup E 



YN,e,n _ Y 



A;=0 '--'tfc 



z 



N.e,n 
tk 



It is easy to see that 



e{N, e, n) ^ C(ei(iV) + e2{N, e) + e-siN, e, n)), 
with ei(A^) and 62 (A^, e) defined in Proposition 14. 2[ and 

2n-l 



e3{N, e, n) := sup E 

0<fcs£2n 



tk tk 



+ j:e 



k=0 



tk 



Zt 



dt 



^N,e,n ^^^^ 



dt 



4.2 Study of the time approximation error 63 (iV, e, n) 

We need an extra assumption. 

(HYl). There exists a positive constant Kf^t such that Vt, t' G [0, T], Vx G M'^, Vy G M, Vz G 

\f{t, X, y, z) - fit', X, y, z) I ^ Kf^t \t-t'\^'\ 

Moreover, we set e = Tn""" and N = n*, with a,b ^ R+'* two parameters. Before giving our error 
estimates, we recall two technical lemmas that we will prove in the appendix. 

Lemma 4.3. For all constant M > there exists a constant C that depends only on T, M and a, such 
that 

2n-l 



JJ (1 + Mhi) ^ C, Vn G N* 



1=0 



Lemma 4.4. For all constants Mi > and M2 > there exists a constant C that depends only on T, 
Ml, M2 and a, such that 



n (l + Mih, + M2-^\ ^ Cn'^''\ 



Firstly, we give a convergence result for the Euler scheme. 

Proposition 4.5. Assume (HXO) holds. Then there exists a constant C that does not depend on n, such 
that 



sup E 

0^fc^2n 



\^tk-xi 



< c- 



n 
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Proof. We just have to copy the classical proof to obtain, thanks to Lemma 1431 

1 2" 



sup E 

0^fes£2n 



^ C sup hi = Chf). 

0s£isS2n-l 



But 



ho = r(i - n-'^/") ^ C 



Inn 



n 



because (1 — n "^Z" ) ~ C'T^^ when n — )• +oo, so the proof is ended. □ 
Now, let us treat the BSDE approximation. In a first time we will study the time approximation error 
on [T-e,T]. 

Proposition 4.6. Assume that (HXO), (HYO) and (HYl) hold. Then there exists a constant C that does 
not depend on n and such that 



sup E 

n</c<2n 



2n-l 



k=n 



■k+i 



tk t 



dt 



Clnn 
^i-2b • 



Proof. The BSDE (14.11) has a linear growth with respect to z on [T — e, T] so we are allowed to apply 
classical results which give us that 



sup E 



tk tk 



2n-l 

+ 5:e 



k=n 



tk+i 



tk 



tk t 



2 




dt 


^ c (e 



n\\2 



^ , Inn 



by using the fact that (7Ar is A^-Lipschitz and by applying Proposition 
Remark 4.7. 



□ 



When a ^ 1 - 26, then e = Tn-" = o{n^^-^ Inn). So we do not need to have a discretization 
grid on [T — e, T].- n + 2 points of discretization are sufficient on [0, T]. 

When a < 1 — 2b, then is is possible to take only [n'^] discretization points on [T — e, T] with 
a + c= 1 — 2b. In this case the error bound becomes 



sup E 



YN,e,n _ yAf,e 
tk ti^ 



2n-l 

+ j;e 

k=n 



tk+i 



tk 



tk t 



to 




dt 





Inn 



+ 



1 



j^l— 2b fi''+'^ 



and the Proposition 14. 61 stay true. 

Now, let us see what happens on [0, T — e]. 

Theorem 4.8. Assume that (HXO), (HYO), (HYl) and (HXl) or (HXl ') hold. Then for all r] > 0, there 
exists a constant C that does not depend on N, e and n, such that 



sup E 

0^fc^2n 



yJV,£,n _ yA^.e 
th tk 



2n-l 

+ ^E 

fc=0 



tk+i 



tk 



tk t 



dt 



C 



l-2b-Ka' 



n 



with K = 4(1 + r])L)^^Ml^. 
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Proof. Firstly, we will study the error on Y. From (14.1b and (14.41) we get 



E 



''tk 



We introduce a parameter 7^ > that will be chosen later. Thanks to Proposition 14.11 and assumption 
(HYO), / is Lipschitz on [tk,tk+i] with a Lipschitz constant := + .-^,2 where = 

2Lf^zMz^2- A combination of Young's inequality + ^ (1 + 7^/1^)0^ + (1 + and properties 

of / gives 



E 



tk tk 



^ (l + 7fc/ifc)IE 



E. 



tk 



Ik 



tfc+i ik+i 
ik+i 



tk 



tk+i 



+Cihk + ^){hl+ l~^'^^E\X,-Xlfds+ I E 



tfe 



f-tk+l 




tk 


'tk 








ds 


■ 


tk + 1 





ds 
(4.5) 



We define 



o rrN,e,n 

So, Z,; 



tk+i\^tk >■ 

Pt^,^ is 1-Lipschitz, we have 



(Zj^'^'"). Moreover, Proposition 14. II implies that Z^'^ = pt^,^{Z^'^), and, since 



^ tk 



Ptk+ii^s ' )-Ptk+iiZt, ) 



<: 



Z 



N,e yN,£,n 



tk 



(4.6) 



As in Theorem [331 we define Z^''^ by 



E,, f Zf'^d. = E,, ({Y.^ll + f"^' fis,X,,Yf'^,Z^'^)ds)\Wt,^, - Wt,, 



tk 



tk 



Clearly, 



E 



tk 



Z 



N,e _ ^N,e,n 



tk 



ds = E 



tk+i 



tk 



^ tk 



ds + /ifcE 



tk *fe 



(4.7) 



The Cauchy-Schwartz inequahty yields 



and consequently 



tk tk 



^ (l + r/)i/3E 
+C/ifcE 



l-k \ IkMl ifco-l 



N,£,ns 



|/(s,x„y/'^zf'-)|'(is. 



(4.8) 
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Plugging (14771) and (l48l) into (1431) . we get: 



E 



N,e 



tk+1 



+ {1 + r])Klihk + —)E 



+C{hk + 



Ik' 



Ik 



j-tk+l 




Itk 





ds 



tk+1 



tk 



'ds + 



■k + l 



E 



tk 



N,e 



Y 



N,e,n 

k+l 



Ik 



^tki 



tk+1 tk+1 



^tkiY,^i[-Y,'''^-- 



ds 

2 



Cfe + 1 



7fc Jtk 



ds. 



Now write 



E 



Y^ 



N,e 



Y, 



N.e,n 



tk+1 



^ 2E 



N,e 



Y 



N,e 



tk+1 



+ 2E 



tk+1 tk+1 



E \Xs - X^X ^ 2E \X, - XtX + 2E jX*, - X 



tk 1 



(4.9) 
(4.10) 



we obtain 

E 



tk tk 



^ (l + 7fc/ifc)lE 



E. 



tk 



Y, 



N,e 



tk+1 
1 



Y 



N,e,n 

k+l 



Ik 



rtk+i 


^ tk 


'tk 





1 



+c(/ifc + -) ( /li + 

Ik 



tk+1 



tk 



+C{hk + 



Ik 



tk+1 



E 



ds 

E \Xs -Xtfds + /ifeE - XI 

2 



7V,e 



Y, 



Ik 



^tki 



tk+1 



ds + /ifcE 

2 



y 



N,e,n 



tk+1 



yiv.fc yA'',e,n. 

tk+1 



N,e,n^ 
tk+1 > 



+CKl{ht, + -)hi,E [ |/(s,X„y/•^Zf'^)| 
Taking 7^ = (1 + r/)^/'^K|: for /i^ small enough, it gives 



ds. 



E 



tk 



Y 



N,e,n 



tk 



^ (1 + Chk + (1 + ??)2/3i^|/ifc)E 



y 



y 



+CE 



+c 



tk 

k+l 



s tk 



tk+1 tk+1 
tk+1 



+ Chi + Chk max E - X 



E|X, 



s - Xt^ I ds 



tk 



E 



y 



N,e 



Y 



N,e 



ds + ChkE 



■k+l 



f{s,Xs,Y,''^',Z^'')^ds 



tk 
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because K'^hk ^ C{ho + hk{T — tfc+i) ^) ^ The Gronwall's lemma gives us 



E 



n-1 
3=0 

+E 



■i-i 



[](l + C/ii + (l + r?)2/3i^2/,. 



i=0 

j+1 



/i^ + /i,- max E X., - 



n|2 



ds 



+ 



[];(l + C/ii + (l + r?)2/3i^2^,) 



E 



Then, we apply Lemma l44l 

2 



E 



ho + maxE Xt, - X?\ 



'j+l 



S tj 



+ x,-Xtr + 



+ ^E 

+/ioE /*" |/(s,X„y,^'%Zf'^)|'(is + E 



Y ''' - Y 



ds 



YN,£ YN,£,n 



A classical estimation gives us E 



Xs — Xt- 



^ \s — tj\. Moreover, since Z^'^ is bounded. 



/o 



^ C7r(l + |y^'"| ) + Cn^''E 



T 



But we have an a priori estimate for E 



Jo 



N,e 



ds 



that does not depend on N and e. So 



E / |/(s,X„y,^'^Zf'^)|'ds^Cn^ 

'0 



2b 



With the same type of argument we also have 

E 



S tj + l 



^ ChjU 



2b 



(4.11) 



(4.12) 



If we add Zhang's path regularity theorem [331 Proposition I4.5l and Proposition 14.61 we finally obtain 

Inn 



E 



n ^^l-2b-{l+^){K^)'^a- 

Now, let us deal with the eiTor on Z. First of all, (14.61 ) gives us 



(4.13) 



n-l 
fc=0 



dt 



71-1 



fc=0 



2N,£,n ^^'^ 
tk t 



dt 
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For ^ /c ^ n — 1, we can use (14.71) and (14.81 ) to obtain 

ik+i 



E 



2 




r-tk+i 


2N,£ ^N,e 

tk t 


2 


dt 


^ E 




dt 






Jtk 





+ (1 + 7?)2/3e 

ftk+i 



tk+1 tfc+1 



)- Ea<;:-<;r) 



+C/ifcE 

Inequality (14.111 ) and estimates for Z give us 



ds 



tk 



n-l 



^N,£,n ^^^^ 
tk t 



dt 



2N,£ _ ^N,e 
tk t 



dt 



A;=0 
n-l 

fe=0 

+(l + r?)2/3^E 

fc=0 

JO 



(4.14) 



+(7/ioE 

n-l 



n-l 



E. 



yN,e _ yN,e,n . _ £ (yA^.^ _ yA^,^>") 



tk+i 



tk 



k=0 

+ (1 + 77)2/3 ^E 



tk t 



dt 



n-l 



A;=0 



E. 



tk^ 



+CE 



tji tji 



*fc tk 



rN,£,n^ 
-tk+i ' 



(4.15) 



with an index change in the penultimate line. Then, by using (14.51 ) we get 



(1 +^)2/3]E 



^t,{ 



Y^ Y-^'^'^ 
tk tf;. 



1 



tk+1 *fc + l 



2 o 1 

+ (l + 77)K|(/ifc + -)E 



7fc Jtfe 



^Af,e _ ^N,£,n 

S tf^ 



+C{hk-\ )hk\hk+ sup E 

\ s€[tk,tk+i] 



YN,£ _ -yA'',e,n 

s tk+1 



ds 



(4.16) 



Thanks to (l4i9l) . (14.101 ). (14.12b and a classical estimation on E 



1^. - 



tk I 



sup E 

se[tfc,tfc+i] 



vr^' - y 



^ C I /ifcn^'' + E 



we have 
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Let us set 7^. = 3(1 + ?7)K|. We recall that /ifciTf < when n 0. So, for n big enough, 

(14.161 ) becomes 



(1 + ??)2/3e 
Clnn 



-E 



n 



y7V,e _ y7V,e,n 



+ -E 
2 



-fc+i 



S th 



ds 



If we inject this last estimate in (14.151 ) and we use Theorem l3.51 we obtain 



k=o '--'tt 



fc+i 



th t 



dt 



^ Chon^^ + C\un sup E 



By using (14.131) and Proposition |4j6j we finally have 



sup E 

0sSfcs;2n 



2n-l 
fc=0 



tk + 1 



tk 



dt 



^ C- 



(Inn) 



i-2b-_fs:a' 



n 



with ii' = 4(1 + 2- Since this estimate is true for every > 0, we have proved the result. □ 

4.3 Study of the global error e{N, e, n) 

Let us study eiTors ei(A^) and 62 (A^, e). 

Proposition 4.9. Let us assume that (HXO) and (HYO) hold. There exists a constant C > such that 

C 



n 



2a-4b ■ 



Proof. We just have to notice that 

is bounded by Cn^. 

For we use the classical Lipschitz approximation 



f{t,Xt,Yf'\Z^'')-f{t,Xt,Yf'\0) 



^ C 



N,e 



and 



Z 



N,e 



□ 



gN{x) = inf ^g{u) + N\x-u\\u£ R'^j . 



Proposition 4.10. We assume that (HXO) holds and g is a-Holder. Then, there exists a constant C such 
that 

ei(iV)^ ^ 



2 ha 

n 1-°^ 



Proof. gpf is a A^-Lipschitz function and gN ^ 9 when N — )• +00 uniformly on M*^. More precisely, 
we have 

c 



N~ 



□ 
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Remark 4.11. For some explicit examples, it is possible to have a better convergence speed. For exam- 
ple, let us take g{x) = (|x|" Ix^o) ^ C and assume that a is invertible. Then, we can use the fact that 
this function is not Lipschitz only in 0, and obtain 



ei(iV) ^ 



C 



2ab 



Xt G 



1/9 



c 



Remark 4.12. It is also possible to obtain convergence speed when g is not a-Holder For example, we 
assume that a is invertible and we set g{x) = HiLi ^xi>o{x). Then 

1/9 



ei(iV) ^ C 



^P((XT).G[0,l/iV]) 



,i=l 



c 

mil 



c 



n 



Now we are able to gather all these errors. 

Theorem 4.13. We assume that (HXO), (HYO), (HYl), and (HXl) or (HXl') hold. We assume also that 
g is a-Holder Then for all rj > 0, there exists a constant C > that does not depend on n such that 

C 



e{n) := e{N,£,n) ^ 



with K = 4(1 + r])Ll^Ml^ 



Proof. Thanks to Theorem l4.8[ Proposition 14.91 and Proposition 14. 1 01 we have 

c c c 



e(n) -J 



-2b-Ka 



+ 



Then we only need to set a := and b := (^2-a){2+K)-2+2a 

Corollary 4.14. We assume that assumptions of Theorem \4.13\ hold. Moreover we assume that f has a 
sub-quadratic growth with respect to z: there exists < /3 < 1 such that, for all t G [0, T], x G W^, 

y £R, z,z' e R^""^, 

\f{t,x,y,z) - f{t,x,y,z')\ ^ {Kj^, + Lj,,(|z|^ + |zf )) \z - z'\ . 

Then we are allowed to take K as small as we want. So, for all r] > 0, there exists a constant C > that 
does not depend on n such that 

, ^ C 

e[n) ^ 



n 

l-a 



2a-46 



+ 



2ab 



to obtain the result. 



□ 



l-a 



Remark 4.15. We are able to specify Remark \477\ in our case, when a = and b = j2 ~a){2+K)~2+2a - 

• When K ^ ^E^- that is to say, when a < 2/3 and K is sufficiently small, then we do not need to 
have a discretization grid on\T — e^T]. 

• When K > V*^" . then it is possible to take only [n"^] discretization points on [T — e, T] with 



1 + 



3a 



{2-a){2 + K) -2 + 2a' 



Theorem l4.13l is not interesting in practice because the speed of convergence depends strongly on K. 



But, we just see that the global enw becomes e(n) ^ ^^1^ when we ai^e allowed to choose K as small 
as we want. Under extra assumption we can show that we are allowed to take the constant M^_2 as small 
as we want. 
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(HX2). b is bounded on [0, T] x M'' by a constant Mf,. 

Theorem 4.16. We assume that (HXO), (HYO), (HYI ), (HX2) and (HXl ) or (HXl ') hold. We assume also 
that g is a-Holder. Then for all r] > 0, there exists a constant C > that does not depend on n such that 

en ^ 



n' 



a—r] 



Remark 4.17. With the assumptions of the previous theorem, it is also possible to have an estimate of 
the global error for examples given in Remarks \4.11\ and \4. 1 21 When g{x) = {\x\'^ l^^o) ^ C, we have 

C 



ein) < 



n 1+29 ' 



and when g{x) = nf=i ^Xi>o{x)> we have 



e(n) < 



C 



Proof. Firstly, we suppose that / is differentiable with respect to z. Thanks to Remark [33] we see that 
it is sufficient to show that 



\Z. 



N,e\ 



ds 



Ft 



is small uniformly in uj, N and e when t is close to T.We will obtain an estimation for this quantity by 
applying the same computation as 121 for the BMO norm estimate of Z page 831. Thus we have 



\z. 



ds 



^(y/"'') +CiT-t), 



with (p{x) = (^Q^<^i.^+^) —2c{x+m) — l) / (2c^), m = \Y\^ and cthat depends on constants in assumption 
(HYO) but does not depend on V^/. Let us notice that m, c and so ip do not depend on N and e. Since 
Y is bounded, 93 is a Lipschitz function, so 



Z. 



ds 



-'T -"t 



Ft 



+ C{T-t). 



We denote by {Y 



N,e,t,x ^N,e,t,x 



the solution of BSDE (14.11 ) when X.'^ = x. As usual, we set X. 



7-t,X 



X 

and z^'^'^'^ = for s ^ t and we define x) := y^^'"^'*'^. Then we give a proposition that we will 

prove in the appendix. 

Proposition 4.18. We assume that (HXO), (HYO), (HYI), (HXl) and (HXl) or (HXl') hold We assume 
also that g is uniformly continuous on W^. Then u^'^ is uniformly continuous on [0, T] x and there 
exists oj a concave modulus of continuity for all functions in |n^'^|A^ G N, e > O}." i.e. u does not 
depend on N and e. 



Then 



Tt 



< E*^ 



N,e\ 



+2\Y 

+2|y^'^ 



00 ^\f^ a{s)dWs\>u 



u'''%T,XT)-u^'%t,X, 
\Xt - Xt 



00 ^\J'^ a{s)dW,\>v 
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21 



with dWs = dWs - VJ%s, Xs, Ys^'', Z^'')ds. But, 



^ h{s,Xs)ds+ r VJ%s,Xs,Yj'^',Z^^')ds+ r a{s)dWs 
t Jt Jt 

T 



^ Mh{T-t) + u + C {l + \zf''\)ds 
^ C{T -t) + v + C{T -tf/'^ (^j^ 



' ds 



1/2 



Since w is concave, we have by Jensen's inequality 



^ u:\C\T-t\+v + C{T-tf/'^^^^ 



T 



\Z^''\'ds 



1/2 







l/2\ 









^ uj\C\T -t\+v + C{T - t)i/2]EQ 

^ ^{c\T-t\ + v + C{T- t)^l^ ■ 

But, oiily depends on constants in assumption (HYO), so it is bounded uniformly in 
and e. Moreover, 



^fa{s)dWs 



is independent of Ft so we have by the Markov inequality 

■T 



^\J^ CT{s)dWs\>V 



Ft 



a{s)dWs 



> V 



CiT-t)^/"^ 



Finally, we have 



E'- 



by setting v = \T - and E«^" 

When / is not differentiable with respect to z but is only locally Lipschitz then we can prove the result 
by a standard approximation. □ 



Ft 



^ uj[C\T -t\^l'^ + v] +C 



XT-t) 



1/2 



oj[C\T -t\^''^ + \T -t^/^] +C\T -t 



,1/4 



Ft~^ — )■ uniformly in oj, N and e when t ^ T. 



5 Some additional results on the time dependent estimate of Z. 
5.1 What happens if a does not depend on time ? 

We have seen that the key point of our approximation results is the time dependent estimate of Z in 
Theorem 13.21 This estimate needs the technical assumption (HXl). When a does not depend on time, 
this assumption becomes 
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(HXl). b is differentiable with respect to x and there exists A G such that for all rj 

|*77cr*cr*V6(s,x)77| s: A|*^7c^|^ (5.1) 
It is possible to obtain some equivalent assumptions. 

Proposition 5.1. Let us assume that (HXO) hold, b is differentiable with respect to x and a does not 
depend on time. Then, following assertions are equivalent: 

(i) Vx G M'^, Vs G [0,r], lmVb{s,x)a C Imcr, 

(ii) Vx G M°', Vs G [0, T], ker *cj C ker ^aVb{s, x), 

(Hi) 3f : [0, T] X ^ M'^^^ bounded such that \fx G M.'^, Vs G [0, T], Vb{s, x)a = af{s, x), 

(iv) there exist A : [0,T] x M'^ ^ W^^'^ and B : [0,T] R^^"^ such that A is differentiable with 
respect to x, V A is bounded and Vx G M"', Vs G [0, T], b{s, x)a = aA{s, x) + B{s), 

(v) there exists A G such that Vry G M'* 

|*?7(T*cj*V6(s, x)r?| ^ A |*?7ap . (5.2) 

Proof. 

(iii) =^ (1), (iii) =^ (ii), (iii) =^ (v) and (iii) 44> (iv): Trivial. 

(i) =^ (iii): There exists a such that aa\ imo- = -^^^l imo-- Since Im V6(s, x)a C Imo" then Vfe(s, x)a = 

aaVb{s, x)a := o"/(s, x). V6 is bounded so / is also bounded. 

(ii) =^ (iii): The proof is the same when we notice that there exists a such that cra\i^j = Id\]\j with 

M e ker a = M"^. 

(v) =^ (ii): Let 5'j'(]R) denote the space of non negative symmetric matrix. Since a^a G 5j"(M), there 
exists s G <S'j'(M) such that aV = s^. It is easy to show that ker s = ker V and ker s*V6(s, x) = 
ker V*V6(s, x). Moreover, V?7 G and pr/c^ = so we are allowed to assume that 

a G 

If a is invertible or o" = 0, the result is proved. Otherwise, let us consider rji and 772 two eigenvec- 
tors of a linked with eigenvalues Ai 7^ and A2 = 0. (v) gives us, for all a G M, 

|*(77i + a772)cr^*V6(s,x)(?7i + 0772)! = | Ai*r/i*V6(s, x)??! + aAf7?i*V6(s, x)772| 

|*(7?i + a?/2)<7^*V6(s, x)(r/i + ar]2)\ is bounded with respect to a, implying *7?i*V6(s, x)r/2 = 0. 

Since a G S'J'(M), we have ker o" © Imo". So, for all 7]2 G ker o", we have shown V6(s, x)r72 G 
(Imcr)-*- = ker (7, that is to say, a^Vb{s, x)r]2 = 0. This last equality allow us to conclude that 
ker cr C ker fT*V6(s,x). □ 

The assertion (v) allows us to reduce assumptions on the regularity of b. Indeed, when A is only a 
Lipschitz function with respect to x with a Lipschitz constant Ka, then it is possible to do a classical 
approximation of ^ by a sequence of differentiable functions {An)nGN such that | V^„(s, x)| ^ Ka- 
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(HXl"). 0- does not depend on time. There exist A : [0, T] x M'' ^ M'^^'^ and B : [0, T] R^""^ 
such that Vx G W^, Vs G [0, T], b{s, x)a = aA{s, x) + B{s). 

Then we obtain a new version of Theorem l3.2l 

Theorem 5.2. Suppose that (HXO), (HYO) and (HXl") hold. Moreover, suppose that g is lower (or 
upper) semi-continuous. Then there exists a version of Z and there exist two constants C, C £ M"*" that 
depend only in T, Mg, Mf, Kf^x< ^f,y> ^f,z and Lf ^ such that, Vt G [0, T[, 



\Zt\ < c + c"(r-t)- 



-1/2 



5.2 Some examples and counterexamples. 

As we already explain in the introduction, Theorem l3.2l is also interesting for BSDEs with a linear growth 
with respect to z. Indeed, such an estimation has been already proved when a is invertible by using the 
Bismut-Elworthy formula, see e.g. ifTOl . while in our case we do not need to assume the invertibility of a. 
Nevertheless, technical assumptions (HXl) or (HXl") appear. Let us see what could happen when these 
assumptions do not hold. The most simple counterexample is given by Zhang in |[20l . Let us consider 
our SDE in dimension 1. We take: 

• T = 2, 

• 6 = 0, 

• VtG [0,2], a{t) = (l-t)lt<i, 

• / = 0, 

• Vx G M, g{x) = arctan (^y^iji 
Proposition 5.3. VM > 0, 30 ^ ti < t2 < I P(|Z|[ti,t2] | ^ M) > 0. 

Proof. Let us denote u{t, x) := Y"/'^ with {Y^''^ , Z^'^) the solution of the BSDE (12.21 ) where X verifies 
Xt = X. For t G [1, 2], obviously x) = g{x) and Z*''^ = 0. For t G [0, 1[, we have 



u{t,x) = E 



g (^x + j\l - s)dw)j 



Vu{t,x) = / g{y) 2 — 6xp — dy. (5.3) 

V27rcr. Jk at \ 2af 



'2i:at JM. \ ^^t 

with (j| = //(I - s)'^ds = \{l - tf. Moreover, 

- / 9{y) —exp — 

By using the substitution y = atz, we get 

Vu{t,0) = J- / g{atz)zexp ( dz 

I I i/A z \ I z^ 



arctan | a/ -^^^ ) zexp ( -— ) dy. 
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By dominated convergence theorem we have 



1 



1/4 



arctan a 



1/4 



zexp 



z 



dy 



I |2-3/4 

\z\ ' exp 



> 0. 



So, 



CJ- 



3/4 (l_t)l/8 



+00. 



Moreover, the equation (15.31) allows us to show that {t,x) i— > Vu{t,x) is continuous on [0, l[xM. So, 
for all M > 0, 3[ti,t2] C [0, 1[ with ti < t2, 3e > 0, such that Vx G [-e,e], Vi G [ti,t2] we have 
|Vti(t,x)o-(t)| ^ M. Since = Z*'^*" = Vu(t, )o-(t), then 



^ e) > 0. 



□ 



To construct this example, the idea is to stop the noise after the time t = 1 to obtain no random 
evolution in the time interval [1,2]. In this example we need to have a function a that depends on time. 
Let us see if it is possible to construct a counterexample when we are not in this situation. If we suppose 
that (HXO) holds, then the assumption (HXl") is always true in dimension 1. So we will take d = 2 and 

• T = 2, 



Vt G [0,2], Vx = ix\x^) G M2, bit,x) = ^^"^^ ) 



yt G [0,2], a{t) 



1 




• / = o 

With these functions, (HXl) does not hold when /i / 0. The solution of the SDE is 



x^ + W} 



x^ + X 



^ / h{s)ds+ I h{s)Wlds. 
Jo Jo 



As we try to obtain an explosion at time t = 1, we will take for h a continuous function such that 

^i[o,if ^ ^ = 0- Moreover, we assume that Vx = (x^,x^) G M?, g{x) = g{x'^) with 



t,x 



^ : M — > M a bounded lower (or upper) semi-continuous function. As before, we denote u{t, x) := 
with (y/'^, Zl''^) the solution of the BSDE (12.21) where X verifies Xt = x. For t G [1,2], obviously 



x) = g{x) = g{x^) and Z^*'^ = 0. For t G [0, 1[, we have 



u{t^ x) 



E 



g \ X + x 



^ [ h{s)ds+ [ h{s)W}ds 
Jt Jt 



Let us denote H a primitive of h. Then, 



/ h{s)Wlds = H{l)Wl - H{t)W} - [ H{s)dW} 
Jt Jt 

= [Hil) - H{t)]W,' + I\h{1) - H{s)]dWl 
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So, 



u{t, x) 



27rat 



g{y) exp 



2al 



with fjf = [H{1) - H{t)Yt + jl[H{l) - H{s)fds and a* = [H{1) - H{t)]. Moreover, 



du 



Zt = Vu{t,Xt)a = -^{t,Xt), 



and 



du 
dx 



j{t,x) 



at 



9{y) 2 



1\2 



{y - - a tx'^) 
2af 



dy. 



So, By using the substitution y — x"^ — atx^ = atz, we get 



du 
dx 



j{t,x) 



^L- I g{atz + x^ + atx^)zeyiY>[-^\dz, 
2-11(71 Jm. \ 2 



and 



du 



/ |z|exp(-— ]dzi^C 



because \at\ ^ at. Finally, we obtain that Z is bounded on the interval [0, 2] even if we it is possible to 

show that ■§^{t, 0) +oo for well-chosen functions g and This example proves that (HXl") is a 
not necessary but sufficient assumption. To be more precise, when a does not depend on time we do not 
succeed to find an example of BSDE such that we have not the estimate 



Zt s: C + C"(T-t) 



-1/2 



So, the assumption (HXl") seems to be artificially restrictive and unnecessary. But this question remains 
open. 



f2n-\ 



6 Appendix 

6. 1 Proof of Lemma 14.31 

We have, 

2n-l 

J] (1 + Mh,) =(11(1+ ^^^') I ( n (1 + ^^^^ 

i=0 \i=0 / \i=n 

Firstly, 

2n-l , TX" 

JJ (1 + M/i,) ^ M + M-j ^C. 

i=n 

Moreover, for ^ ?' ^ n — 1, 



/i, = ti+i -ti = rn-"'/"(l - e-"^) ^ Tn-^'/^'a 



Inn 



n 



'Take for example h{s) = (1 — s)l3<i and g{x) = arctan jV^ )' 
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thanks to the convexity of the exponential function. So 

n— 1 n— 1 



+ Mhi) ^ n ( 1 + MTan""*/"— j 

_ hm / 1 - (1/n") 
^ exp ( MTa— ( — ^ ' ^ 

^ exp I MTa 



\ n \l - {l/n^"-/"^)) 

In n n'^/" \ 



But, 



n n<^/" - 1 / ' 
Inn n'*/" Inn 1 1 



n n"/" -In a— a 
when n +00. Thus, we have shown the result. □ 



6.2 Proof of Lemma 14.41 

Thanks to Lemma 1431 we have 

rn— 1 / 1 I T\/r -L I T\/r_ h. 



n-1 



.■_n \ ~r -^'^-'2t^_+. , , / :.—n 



So, we just have to show that 



But, 



So, 



1 + M2 ^ \ = 1 + - 1). 



n 1 + A^2^3^ = (l + M2K/"-l) 

j=0 ^ 



exp nhi 1 + aM2 hO 



n \ n^ 



exp ( aM2 In n + O 



V ~ \ n 

when n — +00. Thus, we have shown the result. □ 
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6.3 Proof of Proposition 14.181 

We will prove this proposition as the authors of ^ do for their Proposition 4.2. In all the proof we omit 
the superscript N, e for u, Y and Z to be more readable. Let xq, x'q G M.'^ and toj ^ [0; 
argument of symmetry we are allowed to suppose that to ^ We have 

\u{to,Xo) - u{tQ,XQ)\ ^ \u{to,xo) - u{to,XQ)\ + \u{to,Xo) - u{to,XQ)\ . 

Let us begin with the first term. We will use a classical argument of linearization: 



Y 



to,x'o 



ds 



with 



/^(s Xs*^'^" y^*0'^o ^^"'^'o^ _ y^*0'^'o ^*o,a;o^ 



elsewhere, 



/3. 



ye^^ j^to,a^O -yio,a;o ^to,x'o^ _ j^e^g j^to,xo y^°'^0 ^^0,x:'q-^ 





/^(s J^*0'^0 -y^to,a;o ^io,2:o-j 



elsewhere, 



7s 







elsewhere, 



and dWs := dVFs — Jsds. By a BMO argument, there exists a probability Q under which 1^ is a 
Brownian motion. Then we apply a classical transformation to obtain 



Y 



t0,XQ 



T 
to 



ds 



and 



yt(i,Xo Y 



to ,x'o 



+ 



to 



ds 



with uj a modulus of continuity of g that is also a modulus of continuity for g^- We are allowed to 
suppose that w is concav^ so Jensen's inequality gives us 



\u{to,xo) - u{to,XQ)\ ^ C (uj (e' 



yto,XO -y'^0,x'o 



to 



j^to,xo _ j^to,x'g 



ds 



By using the fact that b is bounded we can prove the following proposition exactly as authors of 13 do 
for their Proposition 4.7: 



''There exist two positive constants a and b such that a;(x) ^ ax + b. Then the concave hull of x i— >■ lu{x) V (l^^i (ax + b)) 
is a also a modulus of continuity of g. 
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Proposition 6.1. 3C > that does not depend on N and e such that Vt, t' G [0, T], Vx, x' G M'^, 
Vs G [0, T], 



i^C \x- x'\ + \t-t' 



1/2 



Then, 



\u{to,Xo) - u{to,x'Q)\ ^ C {UJ {\xo - x'q\) + \xo - Xq\) . 

Now we will study the second term: 



\u{t0,XQ) - U{tQ,XQ) 



Firstly, 



to 



yo 



+ 



^t' 



to 



Moreover, as for the first term we have 



/in f^^'^'^ ( ^*0'^0 _ V*0'^0 



'to 



ds 



and 



Finally, 



and 



^to 



^ c[u[\to-t',\'/') + \to-t',\'^'" 



So u is uniformly continuous on [0, T] x M'^ and this function has a modulus of continuity that does not 
depend on N and s. Moreover, we are allowed to suppose that this modulus of continuity is concave. □ 
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